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ABSTRACT 

We present a physically-based absorption-line model for the spectroscopic study of the inter galactic medium 
(IGM) . This model adopts results from Cloudy simulations and theoretical calculations by iGnat & Sternberg! 
(20071) to examine the resulting observational signatures of the absorbing gas with the following ionization 
scenarios: collisional ionization equilibrium (CIE), photoionization equilibrium, hybrid (photo- plus collisional 
ionization), and non-equilibrium cooling. As a demonstration, we apply this model to new observations made 
with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope of the IGM absorbers at z ~ 0.1877 
along the 1ES 1553+113 sight line. We identify Lya, C in, O vi, and N v absorption lines with two distinct 
velocity components (blue at Zb = 0.18757; red at z r = 0.18772) separated by A(cz)/(1 + z) ~ 38 km s . Joint 
analyses of these lines indicate that none of the examined ionization scenarios can be applied with confidence to 
the blue velocity component, although photoionization seems to play a dominant role. For the red component, 
CIE can be ruled out, but pure photoionization and hybrid scenarios (with T < 1 .3 x 10 5 K) are more acceptable. 
The constrained ranges of hydrogen density and metallicity of the absorbing gas are «h = (1.9-2.3) x 10~ 5 cm -3 
and Z = (0.43 -O.67)Z . These constraints indicate O vi and H i ionization fractions, /ovi = 0.10-0.15 and 
/hi = (3.2-5.1) x 10~ 5 , with total hydrogen column density Nh = (0.7- 1.2) x 10 18 cm" 2 . This demonstration 
shows that joint analysis of multiple absorption lines can constrain the ionization state of an absorber, and 
results used to estimate the baryonic matter contained in the absorber. 

Subject headings: Cosmology: observations — intergalactic medium — quasar: absorption lines — ultraviolet: 
general 



1. INTRODUCTION 

The intergalactic medium (IGM) is believed to be the main 
reser voir of bar yons throughout the history of the universe 
(e.g., IShullI 120091) . In the early universe (z > 3) the IGM 
was domin ated by photoioniz ed gas traced by Lya forest 
lines (e.g., iRauch et alJ Il997h . At z<2, the IGM is ex- 
pected to be a mixture of photoionized and collisionally 
ionized gas produced by radiation of active galactic nuclei 
(AGN) and gravitational shock heating during the forma- 
tion of galaxies, groups, and clusters (e.g., Birnboim & Dekel 
2003; Keres et al. 2009). Cosmological hydrodynamic simu- 
lations of large-scale structure formation predict that, at the 
current epoch, >40% of the baryons may exist in the form 
of a warm-hot intergalactic medium (WHIM) at temperatures 
T ~ 10 5 - 10 7 K in a "Cosmic Web" distributed between the 
galax ies (e.g., ICen & Ostriker] 120061: lOppenhei mer & Dave I 
120081) . 

Conducting a complete baryon inventory for the current 
universe is one of the major tasks of modern cosmology. The 
predicted physical conditions in the hottest parts of the WHIM 
may require detection using X-ra y absorption lines (e.g., 
iDave et aTteOOUlCen & F ang 2006). However, owing to lim- 
ited sensitivity and spectral resolution, X-ray observations 
have so far yielded only a few controversial d e tections (e.g.. 
Fane et all 12001 120071: [Nicastro et all l200l IKaastra et alJ 
2006; iBuote et all 120091) . Current X-ray observatories 
such as Chandra and XMM-Newton have spectral resolu- 
tion (A/AA ~ 400) and detection limits (logAfovn£ 15.5; 
lYao et alj|2009l l2010h . Instruments operating in the far- and 
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near-ultraviolet (UV) wavelength bands have superior spec- 
tral resolution, A/AA - 20,000-43,000. With the larger 
intrinsic line strengths (O vi vs. O vn) and the higher spec- 
tral resolution and sensitivity, the UV spectrographs aboard 
the Hubble Space Telescope (HST) - the Space Telescope 
Imaging Spectrograph (STIS) and the Cosmic Origins Spec- 
trograph (COS) - are sensitive to lower column densities, 
logA/ovi ~ 12.7. These instrumen ts can probe the W HIM 
up to redshifts z^O.l (e.g.. iDanforth & Shulll 12001 l200l 

~|2008). 



ISavage et alJ |200H iTripp et all 12008: Thorn & Ched l200i 
The major advances in low-redshift IGM studies in the next 
decade are expected to come from archival data from STIS 
and the Far Ultraviolet Spectroscopic Explorer (FUSE) to- 
gether with new observations from COS. Indeed, the multi- 
phase nature of the IGM has been revealed through detected 
absorption lines of ions at various ionization states (e.g., H i, 
C m, C iv, Si m, Si iv, N v, O vi, and Ne vm) by FUSE, STIS, 
and COS. In order to proceed, we need to know how to convert 
measurement of these absorption lines to amounts of baryonic 
matter, which depends on ioni zation fractions, metallicity, and 
ionization phases of the IGM dPenton et alj|2004t ITripp et al] 
l2008l:lDanforth & Shulll200l hereafter DS08). None of these 
properties is well measured for many absorbers, and in many 
cases the data do not permit such determinations. 

Spectroscopy of multiple absorption lines provides a pow- 
erful diagnostic of IGM properties. These diagnostics have 
been applied to many systems by measuring ionic col- 
umn densities and doppler velocities, examining the signifi- 
cance of photo-, collisional, and non-equilibrium ionization, 
and assessing the p h ysical conditions in t he absorbers (e.g. , 
[ Richter et all 120041; [Savage et all 120051; ITripp et all 120081; 
Narayanan et al.ll2010l) . However, because these applications 
follow a multistep procedure, the derived physical parame- 
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ters could depend in complex ways on the initially measured 
quantities (e.g., temperature and metallicity of an interven- 
ing gas could co-vary as different functions of O vi column 
density), which makes estimations of their uncertainty hard to 
implement. We present an absorption-line model in which the 
physical properties of the absorber have been used as the pri- 
mary fitting parameters. Spectral fitting of multiple absorp- 
tion lines with this model will not only yield constraints to 
these parameters but automatically take care of covariance of 
all parameters. As a case study, we apply this model to an 
interesting and well-defined absorption system at z ~ 0.1877 
along the sight line of 1ES 1553+1 13 observed with COS. 

The paper is organized as follows. We describe our COS 
observations and data reduction procedure in § |2] and intro- 
duce our spectral model in § [3] We apply our model to the 
observations and discuss our results and their implications in 

m 

2. OBSERVATIONS AND DATA REDUCTION 

One of the Guaranteed Time Observation (GTO) targets, 
the BL Lacertae Object 1ES 1553+1 13 was observed on 2009 
September 22 with COS, which was installed on the HSTdur- 
ing Hubble Servicing Mission 4. The observations were taken 
with the same FP_POS position (FP_POS=3) but with differ- 
ent central wavelengths, shifting along the sequence 1291 A, 
1 300 A, 1 309 A, and 1 3 1 8 A for the G130M grating (covering 
1 134 A to 1465 A), and 1589 A, 1600 A, 161 1 A, and 1623 A 
for the G160M grating (covering 1383 A to 1797 A). The to- 
tal exposure times are 3109.6 s with the G130M and 3802.8 s 
with the G160M. 

The observations were calibrated with the pipeline CAL- 
COS (version 2. 1 If). Flat-fielding, alignment, and coaddition 
of the processed exposures were carried out using IDL rou- 
tines developed by the COS GTO team specifically for COS 
FUV datcQ. Briefly, true fiat-fielding of COS observations has 
yet to be developed, but each exposure was corrected for nar- 
row, ~ 15%-opaque shadows from repellor grid wires. The 
one-dimensional map of grid-wire opacity was shifted from 
detector coordinates to wavelength space and divided from 
the reduced one-dimensional flux. Within the same grating 
(G130M or G160M), these spectra were then cross-correlated 
and combi ned to form an expo sure-weighted coadded spec- 
trum (See iDanforth et al.ll2010bl for a detailed description of 
observations and data reduction). 

In the COS spectra of 1ES 1553+1 13, three absorption sys- 
tems w ere identified near redshift z = 0.187 (Danforth et al. 
2010b). Among these systems, the one at z ~ 0.1877, as pre- 
sented in Figure [T] has the most significant detections of ab- 
sorption lines of Lya, C in, and the O vi and N v doublets. We 
therefore use this system to test our absorption-line model. 



3. MODELS 

We considered intergalactic absorbers in the following ion- 
ization model conditions: (1) collisional ionization equilib- 
rium (CIE; Case 1); (2) photoionization equilibrium (PE; 
Case 2); (3) medium with a constant electron temperature illu- 
minated by the UV photoionizing background (hybrid model; 
Case 3); and (4) post-shock gas undergoing non-equilibrium 

2 See http : / /casa . Colorado . edu/~danf or th/ cos too Is . html 
for our coaddition and flat-fielding algorithm and additional discussion. 



TABLE 1 
Ionization Scenarios 



Scenarios Free Parameters 

Case 1 : collisional ionization equilibrium (CIE) T, Z, Nn 

Case 2 : pure photoionization equilibrium (PE) nu , Z, Nu 

Case 3 : photo- plus collisional ionization (hybrid) njj, T, Z, Nu 

Case 4 : non-equilibrium isobaric cooling T, Z, Nu 

Case 5 : non-equilibrium isochoric cooling T, Z, Nu 

NO TE. — Cases 1, 4, and 5 are based on results by Gnat & Sternberg 
12007), and Cases 2 and 3 are based on Cloudy simulations. «h, T, 
Z, and Nu are the number density, temperature, metallicity, and total 
column density of the absorbing gas, respectively. 

isobaric (Case 4) or isochoric (Case 5) cooling. These five 
cases represent the most explored ionization scenarios in ul- 
traviolet spectroscopic studies of the IGM. 

For Cases 1, 4, and 5, we adopted the results of 

iGnat & Sternberg! d2007), in which they calculated the ionic 
abundances for gas in CIE state and for a shock-heated gas 
cooling from \ogT = 6.7 either at constant density (isochori- 
cally) or at constant pressure (isobarically). 

For Case 3, we ran Cloudy (version 08.00; iFerland et al.1 
1998) simulations with parameters of temperature (T) and 
hydrogen numbe r density (»h) - A 100-kpc thick slab (e.g., 
iDave et alj 1200 U iPenton et ail 120041) was used to approxi- 
mate an intergalactic absorbing clo ud, which is illum i nated by 
the intergalactic ionization field of Haardt & Madaul (120011) at 
redshifts z = — 1 with a step size of 0.1. In our spectral 
fitting described below, the estimated ionization field at the 
closest redshift bin to that of the observed absorber was uti- 
lized. The grids of log T and log[«H(cm" 3 )] were 4 to 7.5 and 
—6 to -2, respectively, with a step size of 0. 1 for each grid. 
The collisional ionization effects were included by setting a 
constant electron temperature in the simulations. The rela- 
tive ionic column densities were initially calculated with the 
metallicity log(Z/Z Q ) = -1, and the results for metal (heavier 
than helium) ions were then scaled to different metallicities. 
We took the results with log T = 4 for a pure photoionization 
model (Case 2). 

Table [T] summarizes these ionization scenarios, and Fig- 
ures |2] and [3] demonstrate the differences in O vi and N v 
ionization fractions among different models. In any of these 
models, for a given set of physical parameters (e.g., Z, T, 
and/or «hX there will be a determined set of column density 
ratios among all ions. In turn, if the observed multiple absorp- 
tion lines occur in the same gas undergoing the same ioniza- 
tion scenario, the column density ratios among metal ions will 
provide constraints on T and/or riu, and ratios between metal 
ions and H i will yield metallicity diagnostics of the absorbing 
gas. 

We developed a model that can be used to jointly ana- 
lyze multiple absorption lines to realize these diagnostics. 
The model, which is revised from ou r X-ray absorption- 
line model absline (lYao & Wangl 120051) . is summarized as 
follows. An absorption-l ine profile can be described as 
dRvbicki & Lightman|[l979l) : 

/(A) = / c (A)exp[-r(A)], (1) 

where f c is continuum flux and r is the absorption opti- 

3 We in fact used the command "HM05" in Cloudy simulations to spec- 
ify both shape and flux of the ionization field; the normalization of the field 
evolves from 1.6 X 10~ 6 erg s -1 cm -2 at z = to 1.3 X 10~ 5 erg s -1 cm -2 at 
Z=l. 
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FIG. 1 . — Absorption lines at z ~ 0.1877 detected in the COS spectrum of 1ES 1553+1 13. The vertical eiTor bar shown at —270 km s~' in each panel indicates 
the averaged flux error in the plotted wavelength range. Thick black lines show the continua and total absorption, and blue and red lines mark fits to the two 
velocity components, at zb = 0.18757 and z r = 0.18772. The labeled line shifts (v c ), doppler velocities (&d), and column densities (JVi) are measured relative to 
z = 0. 1 8757. Numbers in brackets "[]" are calculated values with respect to other ions. Doublet lines are jointly measured when applicable. Data are binned by a 
factor of 2 with ~ 3 bins in each resolution element. See text for details. 
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FIG. 2. — Ionization fraction of O VI under different physical conditions. 
The solid-blue, dashed-red, dash-dotted-green, and dark-red-dotted curves 
indicate the CIE, Isobaric, Isochoric, and Hybrid models, respectively. For 
the hybrid model, plots show the fraction in gas illuminated by the photoion- 
ization field at z = 0. 
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FIG. 3. — Same as Figure|2]but for ion N V. 

cal depth, given by a Voigt function depending on wave- 
length (A), column density (TV;), doppler width (bo), oscillator 
strength (f u ; / and u are the lower and upper excitation levels 
for the transition), and the natural damping factor (7) of the 
transition. For a set of input parameters (A_u, Z, T, and/or n H ) 
of any ionization case discussed above, the model queries the 
simulation and/or calculation results and interpolates the val- 
ues lying between the simulated and calculated grids to obtain 
a unique value for Nj. These physical properties are fitting pa- 
rameters linked together in a joint analysis of multiple lines. 

We implement our spectral analysis within the software 
package Xspec0(Arnaud 1996), to utilize its % 2 minimization 
algorithm over multiple fitting variables, where 




Here, for each wavelength bin (;'), ff, ff, and c, are the ob- 
served flux, the model-predicted flux, and the uncertainty of 

4 http://heasarc.nasa.gov/docs/xanadu/xspec/ 



ff. For COS gratings G130M and G160M, we interpolate the 
line spread functions^ (LSFs) to the desired wavelength grid 
and then use the IDL scripts wrt_ogip_rm^ and writeph£\ 
to transform these LSFs and the coadded spectra obtained in 
§ [2] to a "fits" format for input to Xspec. In spectral fitting, 
the spectral models are automatically convolved with these 
wavelength-dependent LSFs. 

In the following, we apply this model to the absorption 
system at z = 0.187 along the sight line to 1ES 1553+113. 
Because velocity components in this system can clearly be 
resolved (Fig. [T}, the spectral analyses conducted here are 
similar to traditionally employed profile-fitting techniques, al- 
though our physical variables are the fitting parameters. When 
several velocity components are severely blended, arbitrarily 
resolving components through spectral fitting could lead to 
unreliable results. 

4. SPECTRAL FIT, RESULTS, AND DISCUSSION 

Using our absorption-line model, we first fit the lines indi- 
vidually to obtain properties of the absorbing ions. The local 
continua are interpolated from the best fits to spectral ranges 
±150 km s" 1 from line centroids. Absorption lines of C in, 
O vi, and N v clearly show two velocity components (blue 
and red) in the z ~ 0.1877 absorber, but component separa- 
tion is not obvious in the broad Lya, other than asymmetry 
in the red wing (Figure [TJ. In both components, centroids of 
the metal lines are well aligned, except for the red N v that 
deviates ~ 20 km s" 1 from others. We therefore set the cen- 
troid velocities of the two H 1 components equal to those of 
the corresponding C in absorption. Fitting these lines individ- 
ually yields the velocity shift v c , bo, and N{ of the absorbing 
ions. For each velocity component, we link bo, N, and the 
centroids of the O vi and N v doublets to get better constraints 
(except the blue component of O vi). The O vi 1037.62 A 
line should be about half the strength of O vi 1031.93 A, but 
it appears stronger in the blue component. This mismatch is 
likely caused by an unidentified weak intergalactic Lya ab- 
sorber at z = 0.01364 that contaminates the blue component 
of O vi 1037.62 A. We therefore only link the bo and veloc- 
ity centroids and vary Afovi independently in the doublet. We 
also exclude this line in our joint analysis of multiple absorp- 
tion lines described below. The red wing of the 1031.93 A 
line is blended with Ly/3 absorption at z = 0.1950. We use a 
two-velocity model to fit the Ly/3 profile and fix its contribu- 
tion in our analysis (Fig. [TJ. The bo of C 111 and N v cannot be 
constrained in the blue component, owing to the narrowness 
of the C m line and the poor S/N around the N v line. We set 
the bo ranges for C in and N v to those of H 1 and O vi as- 
suming thermal broadening. Because there is no visible Si m 
1206 A line in the spectrum (see iDanforth et aill2010bl) . we 
fix the line centroid velocities to those of C m and bo to the 
corresponding thermal values of H 1 to obtain upper limits on 
Nsm- The constrained values of v c , bo, and Ni are labeled in 
Figure Q] 

Next, we jointly analyze these absorption lines in each ve- 
locity component to examine the five ionization scenarios dis- 
cussed in § [3] Line widths of IGM absorbers have been 
shown to possess significant contributions from turbulent 
motion or unresolved components; measured widths there- 

5 http://www.stsci.edu/hst/cos/performance/spectral_resolution/ 

6 http://hea-www.harvard.edu/PINTofALE/pro/util/ 

7 http://astro.uni-tuebingen.de/software/idl/aitlib/fits/ 
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fore pr ovide upper limits on thermal broadening in the ab- 
sorber (iThom & Chenll2008tlTripp et al.ll2008l iDanforth et al.l 
l2010al) . Attributing the constrained bo to purely thermal 
broadening, we derive upper limits of log T < 4.7 and 5.1 for 
the absorbing gas in Cases 1 and 3-5 (both are based on H i 
measurements) for blue and red components respectively. We 
apply the upper limits in our joint analyses, and constrain the 
physical properties (Z, T, «h) and total hydrogen column den- 
sity A^h of the absorbing gas (Tabled. Given the constrained 
values of Z, T, «h, the column density of any ion can be calcu- 
lated. The predicted column densities log Nsuu = 8.82- 1 1 . 16 
for all acceptable scenarios, co nsistent with the upper limits 
derived from the data (Figure [T] Danforth et al. 2010b). 

The spectral fitting indicates that some ionization scenarios 
can be ruled out by the data. To quantify how well these sce- 
narios fit the velocity components, we record the difference 
(6x 2 ) °f their best-fit \ 2 from those of fitting these lines indi- 
vidually. According to an F-test, values S\ 2 < 0.1, 1.0, and 
4.0 indicate that a joint fit is identical to an individual-line fit 
at > 94% (99%), 55% (75%), and 9% (19%) confidence levels 
for Case 3 (Cases 1, 2, 4, 5), respectively. In this work, a joint 
fit is considered acceptable if it is identical to the individual- 
line fit at > 10% confidence level. None of the examined 
ionization scenarios is acceptable for the blue component, al- 
though photoionization seems to play an important role in the 
absorbing gas. The non-fruitful exercise for the blue com- 
ponent suggests that the low and high ions are not produced 
in the same absorber in the same ionization conditions (see 
discussion below). As an alternative, non-equilibrium cool- 
ing plus photoionization could reconcile all the measurements 
(J. Li 2010, private comm unication). Unfortunately, this sce- 
nario was not calculated in lGnat & Stern berg (2007). For the 
red component, PE and hybrid ionization cases (Cases 2 and 
3) are acceptable at >75% confidence levels (Table |2), the 
non-equilibrium cooling scenarios (Cases 4 and 5) are less 
acceptable, and the CIE scenario (Case 1 ) can be ruled out. 

Besides the ion column density ratios, differences in ther- 
mal broadening of absorption lines of different elements also 
constrain the gas temperature. We assume that the line 
broadening consists of both thermal and non-thermal (tur- 
bulent) broadening, related by b 2 D = (bfo er + b' 1 mb ). Taking 
kT = 1 /2mb 2 her , we can use the measured bo in different ele- 
ments to determine T, b t h eT , and b tur b- For the red component, 
measurement of H i and C in (Fig.[TJ yields 4.9 < logT < 5.1 
and km s" 1 < b tm \, < 16 km s -1 , while the combination of 
H i and O vi yields 4.8 < log T < 5 . 1 and 1 1 km s" 1 < b tm b < 
22 kms" 1 . These constraints (particularly on temperatures) 
are consistent with those obtained from the column density 
ratios (Tabled. 

In our analysis, we assumed that both low and high ions 
were produced in a co-spatial absorber with the same ioniza- 
tion conditions. This assumption is often justified by the simi- 
larity in their absorption-line kine matics (Figure [p. How ever, 
it has also been suggested (DS08; lNaravanan et al.ll2010l) that 
these ions may be multi-phase (non co-spatial), with the low 
ions (e.g., H i, C in, Si m) arising from cool-phase photoion- 
ized gas and the high ions (e.g., O vi and N v) arising pri- 
marily from collisional ionization. Such scenarios could oc- 
cur, for instance, in IGM filaments with p hotoionized cores 
surrou nded by infalling shock-heated gas (Danforth & Shulll 
200$. To examine this possibility, we use the PE model (Case 
2) to jointly fit the blue Hi, C m, and Si m lines where all our 
models failed in joint analyses of high and low ion lines (Ta- 
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FIG. 4. — The lcr co-varying contours of log (Z/Zq) vs. log tin for pho- 
toionization model, obtained by jointly fitting the blue H I, C III, and Si III 
lines and using the observed amount of O VI (dotted) and N V (dashed) as 
upper limits. 

ble|2]). Because we have only two detections (H i and C m) and 
one upper limit (Si m) and there are three free parameters («h, 
Z, and Nh', Tabled, we cannot obtain meaningful constraints 
on these parameters. Requiring the amount of O vi and N v 
produced in the PE ionization to be less than those observed, 
log AWi < 13.52 and log AW < 13.30 (Fig.Q), we obtain the 
lcr co- varying contour of log (Z/Z Q ) vs. log n# (Fig. The 
overlapped region in Figure|4]indicates the constrained ranges 
of log n H = (-4.65,-3.87) and log (Z/Z Q ) = (-0.65,0.25) for 
the cooler PE ionization gas. While the amount of O vi and 
N v contained in the PE ionized gas is still poorly constrained, 
it is hard to use the remaining amount to further constrain the 
properties of the hotter collisional gas. Once multiple ioniza- 
tion states of low and high ions are available, we can use this 
method to constrain properties of different phases and probe 
the interactions between them. 

All acceptable scenarios require the metallicities of the ab- 
sorbing gas to be >40% of the solar value (Table O. Such 
high metallicities suggest that we are sampling enriched in- 
tragroup gas or a circumgalact ic medium with a small impact 
distance to a galaxy (e .g., iJenkins et al.l 120051; lAracil et"afl 
2006; Prat t et all 120071) . A deep galaxy survey along the 
1ES 1553+113 sight line (B. Keeney et al. 2011, in prepa- 
ration) will examine these possibilities and suggest a physical 
origin for these absorbers. As the metallicity increases, the 
cooling times will become shorter compared to the recombi- 
nation time. The absorbing gas could then become an "over- 
cooled" plasma in which the high ions (C iv, N v, O vi) are no 
longer good ind icators of the gas tempe rature, especially in 
cooling models (Gnat & Sternberg 2007). In this case, tem- 
perature diagnostics can still be obtained from the differential 
temperature dependence between low ions (e.g., H i and C m) 
and high ions, if they are co-spatial, or from differential ther- 
mal broadening of different elements. 

We assum e d rel ative solar metal abundances of 
lAsphind et all d2009h in our spectral analysis. This as- 
sumption can be relaxed when absorption lines from different 
ionization states of the same elements are available to con- 
strain the physical parameters of the absorber. Spectroscopy 
of multiple absorption lines from different elements will then 
yield their relative abundances. 

The joint analyses presented in this work also enable us to 
directly estimate the total baryon content of Lya and metal- 
line absorbers and the uncertainty in those measurements. 
The primary systematic uncertainties arise from modeling the 
metallicity and ionization fractions. Since our model uses 



6 



Yao et al. 



TABLE 2 
Diagnostic Results 





S X 2 


Blue component zt 
log« H log 7" 
(cm" 3 ) (K) 


= 0.18757 

log(Z/Z ) 


log^H 

(cm" 2 ) 


s x 2 


Red component z r — 
logn H log T 
(cm" 3 ) (K) 


0.18772 
log(Z/Z ) 


log^H 

(cm" 2 ) 


Case 1 
Case 2 
Case 3 
Case 4 
Case 5 


MOO 
21.5 
16.6 
MOO 
MOO 


NA 

_ 4 6 n+0.04 

-4 62 mos 
H - -0.08 

NA 
NA 


[<4.7] 

NA 
[<4.7] 
[<4.7] 
[<4.7] 


X 

05+* 09 

u - UJ -0.08 

07+0 07 

X 
X 


X 

17 51+* 05 
17 57« 09 

1 '- J '-0.08 
X 

X 


>100 
0.8 
0.1 
4.1 
3.9 


NA 
_4 74+0.03 

"-0.03 
_4 67+0.03 

NA 
NA 


[<5.1] 
NA 

[<5.1] 
(>5.0), [<5.1] 
(>5.0), [<5.1] 


X 

-0 22+* 05 
u - zz -0.05 

-0 24 m07 

o n+o '4 

a Q9+0.06 


17.89+°;° 3 
17.96* 
18.58«; 
18.59+™ 2 



NOTE. — Cases 1-5 are for gas in pure CIE, pure photoionization, photoionization plus collisional ionization, non-equilibrium 
isobaric cooling, and non-equilibrium isochoric cooling states, respectively. Sx 2 is the difference of the best-fit x 2 of a joint analysis 
from that of fitting lines individually. The "NA" symbol indicates that the parameter is not applicable to the model, and the "x" symbol 
indicates that the model can be ruled out by the data. Limits within parentheses "()" are obtained in joint analysis, and limits within 
brackets "[]" are derived in fitting individual lines. Errors are reported at la uncertainty. See text for details. 



results of simulations and calculations, the metallicity is a 
known parameter and the ionization corrections can easily 
be calculated. Once the physical parameters are constrained 
from spectral fitting, we can infer the column density of to- 
tal hydrogen. For example, the best acceptable photoioniza- 
tion and hybrid ionization scenarios for the red component 
indicate that the absorber contains a total YVh = (0.7- 1.2) x 
10 18 cm" 2 (Table|2]). We are applying this joint analysis tech- 
nique to systematically study intergalactic absorption lines, 
constrain the physical properties of the absorbers, and esti- 
mate the total baryons contained within the absorbers. 

In summary, we have presented a physically-based ab- 
sorption line model for the spectroscopic study of the IGM. 
In this model, we can examine the photoionization equilib- 
rium, collisional ionization equilibrium, photoionization plus 
collisional ionization, and post-shock gas undergoing non- 
equilibrium isobaric and isochoric cooling. For the WHIM, 
in which 40-50% of the baryons at the current epoch are 
believed to reside, non-equilibrium cooling plus photoioniza- 



tion should be a more proper scenario; however, such calcu- 
lations are not yet available in the literature. For the purpose 
of demonstration, we only simulated a 100-kpc thick (e.g., 
iPenton et alj 12002b intervening cloud for the absorbing gas 
along the sight line to 1ES 1553+1 13. Because our model di- 
rectly uses the results from calculations and simulations, any 
progress in future in modeling the WHIM and in simulating 
the 10-100 kpc scales of the circumgalactic medium can eas- 
ily be adopted. 
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